Critical velocity in cylindrical Bose-Einstein condensates 
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We describe a dramatic decrease of the critical velocity in elongated cylindrical Bose-Einstein 
condensates which originates from the non-uniform character of the radial density profile. We 
discuss this mechanism with respect to recent measurements at MIT. 



Superfluidity is one of the striking manifestations 
of quantum statistics, and the occurence of this phe- 
nomenon depends on the excitation spectrum of the 
quantum liquid. The key physical quantity is the crit- 
ical velocity v c , that is, the maximum velocity at which 
the flow of the liquid is still non-dissipative (superfluid). 
The well-known Landau criterion |Q gives the critical ve- 
locity as the minimum ratio of energy to momentum in 
the excitation spectrum: 



e(k) 
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(we put h = 1 and the particle mass M = 1). Actually, in 
liquid 4 He the Landau criterion overestimates the critical 
velocity. The explanation of this fact was put forward by 
Feynman (see e.g. [|J) who suggested that superfluidity is 
destroyed by spontaneous creation of complex excitations 
(vortex lines, vortex rings, etc.). Extensive theoretical 
and experimental studies on this subject are reviewed in 

Bose-Einstein condensation of dilute trapped clouds of 
alkali atoms j|-^3| offers new possibilities for the investi- 
gation of superfluidity . In the spatially homogeneous 
case, the spectrum of elementary excitations of a Bosc- 
condensed gas is given by the Bogolyubov dispersion law 
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and the Landau critical velocity is equal to the speed of 
sound c s = ^JgriQ (no is the condensate density, g — Ana, 
and a > is the scattering length). The first exper- 
imental observation of the critical velocity in trapped 
gaseous condensates, recently reported by the MIT group 
M, gives a significantly smaller value of v c . 

The analyses i n Pi and in recent theoretical publica- 
tions (see e.g. |Tc|jl l[] and references therein) employ the 
Feynman hypothesis and provide a qualitative explana- 
tion of the MIT experimental result. In this Letter we 
point out a simple geometrical effect which is charac- 
teristic for elongated cylindrical traps. Due to the non- 
uniform character of the radial density profile, the spec- 
trum of axially propagating excitations in these traps is 
very different from the Bogolyubov dispersion law (|^), 
and this difference leads to a strong decrease of the criti- 
cal velocity. We show that this effect can at least partially 



explain the small critical velocity measured in the MIT 
experiment 

We consider an infinitely long cylindrical condensate 
which is harmonically trapped in the radial (p) direction. 
Then the condensate wave function V'o(p) satisfies the 
Gross-Pitaevskii equation 



g\Mp)\ 2 MP) = 0, 



where p is the chemical potential, V(p) = u> p /2 is the 
trapping potential, and to the trap frequency. In the 
Thomas- Fermi regime, where the ratio ij = /i/w»l, the 
density profile is given by no(p) = |V>o| 2 = (a* — V{p))l 9 
and the chemical potential is related to the maximum 
condensate density as p — no ma xg- 

Elementary excitations can be regarded as quantized 
fluctuations of the condensate wavefunction ffl. In our 
trapping geometry they are characterized by the axial 
(z) wavevector k and radial angular momentum m. The 
corresponding part of the field operator reads 
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where b m ] i (l)' mk ) are annihilation(creation) operators of 
the excitations. The excitation wave functions can be 
written in the form 

(u, v) mk {p, z) = (ft, v) mk (p) exp(im</>) exp(ifcz), 

where cf> is the angle in the x, y plane, and the radial func- 
tions u m k(p) and v m k{p) are solutions of the Bogolyubov- 
de Gennes equations (see e.g. Jf2|) 



(3) 
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< r = | — ^ + — + V-P + 2n Q gj v + n gu. (4) 

Equations (||) and (||) constitute an eigenvalue prob- 
lem. For given m and fc, they lead to a set of frequencies 
£nm(k) characterized by the radial quantum number n 
which takes integer values from zero to infinity. In the 
limit e nm (k) <C p, these modes were found for m = 



in the hydrodynamic approach in [13-M|. For kR ^ 1 , 
where R = (2/x/w 2 ) 1 / 2 is the Thomas- Fermi radial size of 
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the condensate, the dispersion relation can be expanded 
in powers of k 2 

e 2 n0 (k) = 2cj 2 n(n + 1) + ^(kR) 2 + 0(fc 4 ). (5) 

The lowest mode (n = 0) represents axially propagating 
phonons: For this mode we have eoo(fc) = czk, where the 
sound velocity cz = y riQ max g/2 is smaller by a factor 
of \[2 than the Bogolyubov speed of sound at maximum 
condensate density, c s 13 - 15 1 . The velocity c s /a/2 of 
axially propagating phonons has been measured in the 
MIT experiment @. 

The first correction to the linear behavior of the dis- 
persion law eoo(fc) a t kR <C 1 reveals its negative curva- 
ture: 5e 0Q = -uj(kR) 3 /192 ||. In || the perturbative 
analysis leading to Eq.(|^) was extended numerically to 
k ~ l/R, still assuming that e n o(k) <C The calcu- 
lations show that the group velocity of the first mode, 
deoo/dk, monotonously decreases with increasing k and 
can become significantly smaller than cz characteristic 
for k <C l/R. This indicates that the critical velocity (|l|) 
associated with creating axial excitations (m = 0, n = 0) 
is smaller than c s and can also be reduced to below cz- 
The physical reason is that the decrease of the conden- 
sate density with increasing p makes the axial superfluid 
flow less stable (see below). 

The hydrodynamic approach used in Jl3|— p~5[| is not 
valid for e > /i, where the excitation spectrum is no 
longer phonon-like and is dominated by the single par- 
ticle dispersion relation e(k) = k 2 /2 (see Eq.(|^)). The 
crossover between the two regimes occurs at k k c — 
p}l 2 and prevents the decrease of the group velocity with 
further increase in k. Obviously, the decrease of the criti- 
cal velocity due to the radial inhomogeneity of the density 
profile can be dramatic only for k c 3> l/R, i.e. in large 
condensates with 77 ^> 1. 
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FIG. 1. The dispersion law for the lowest mode (n = 0, 
m — 0) at various values of r\. The excitation energy eoo is 
given in units of /i, and the momentum k in units of k c . 

The Bogolyubov-de Gennes equations (||) and (^) are 



valid for arbitrary k and hence allow us to find the exci- 
tation spectrum in the crossover regime (see Figjl]) and 
establish the value of the critical velocity as a function of 
the Thomas-Fermi parameter rj. Since the spectrum of 
excitations consists of a number of independent branches 
characterized by the quantum numbers n and m, Eq.(Q) 
gives a value i4" m ' for the critical velocity correspond- 
ing to each mode. The results of our numerical calcu- 
lations for the two lowest modes (n = 0, m = and 
n = 1, m = 0) are presented in Fig|2. The breakdown of 
superfluidity occurs when the velocity of the flow matches 
the lowest of the velocities Vc ™* , which is proved to be 
The corresponding curve in Fig.^| indicates a sig- 
nificant decrease of the critical velocity compared to cz 
already at 77 ~ 10. 




FIG. 2. The critical velocity u'"" 1 ' (in units of c s ) for 
the two lowest modes versus the Thomas-Fermi parameter 
•q — /j,/lu. The dashed line shows the value of cz- 

The decrease of the critical velocity with increasing 
ratio fi/co seems counterintuitive. One can increase this 
ratio by decreasing u> and keeping constant fi. Then, at 
a constant density on the axis of the cylinder, for radi- 
ally larger condensates one gets smaller v c - However, this 
phenomenon has a clear physical explanation. The key 
point is the non-uniform character of the radial density 
profile. With increasing axial wavevector k, the wave 
functions of the excitations with k < k c are more local- 
ized in the outer spatial region of the condensate and are 
thus more sensitive to the small value of no in this re- 
gion. It is this feature that provides a decrease of v c with 
increasing 77, since for larger rj one has more possibilities 
to increase k and still satisfy the condition k < k c . The 
described situation is not met in liquid helium where the 
density is practically constant and changes only in the 
region very close to the border of the sample. 

We now turn to the discussion of the MIT experi- 
ment ||. The radial frequency in this experiment was 
u) = 2ir x 65 s _1 and the chemical potential /x = 110 nK, 
so that 77 « 35. The sound velocity c s — 6.2 mm/s was 
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measured by observing a ballistic expansion. For these 
parameters we find v c w 0.42c s = 2.6 mm/s, which is not 
far from the observed value vi expS> w 1.6 mm/s p. We do 
not pretend to explain the MIT data as the condensate 
in the experiment is not an infinite cylinder. The ratio 
of the radial to axial frequency for the MIT trap is 3.3, 
which indicates that the effect of a finite axial size can 
be important. Nevertheless, our calculations may par- 
tially explain the small critical velocity observed in the 
experiment. Moreover, a smaller value of v c originating 
from our geometrical effect can facilitate the nucleation 
of vortex objects and thus provide a larger decrease of 
the critical velocity. 
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